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INTRODUCTION 

'^In  the  past  few  years,  there  has  been  a  renewed  interest  in  the  structure  of  an 
electrolyte  at  an  electrode  interface  (the  so-called  electric  double  layer).  If  only  the  ions  are 
taken  into  account  explicitly,  with  the  solvent  represented  by  a  uniform  dielectric  medium  of 
constant  dielectric  constant  e,  the  potential  difference  across  the  double  layer  is  much  smaller 
than  is  observed  experimentally.  This  remains  true  even  if  the  Poisson-Boltzmann  theory  of 
Gouy^  and  Chapman^  (GC)  is  replaced  by  more  refined  theories,  such  as  the  hypemetted 
chain  approximation.^'^ 

Presumably,  this  error  is  due  to  the  simplified  treatment  of  the  solvent.  If  it  is  assumed 
that  the  dielectric  constant  e  in  the  so-called  inner  region  between  the  hard,  flat,  charged 
wall  (the  electrode)  and  a  plane  passing  through  the  centers  of  the  charged  hard  spheres  in 
contact  with  the  wall  is  different  from  the  bulk  value  e,  ^he  potential  drop  across  the 
interface  is 


V  -  -H  +  ^0'  O) 

where  a  is  the  diameter  of  the  charged  hard  spheres  (o/Z  is  then  the  distance  of  closest 
approach  to  the  wall),  —  E/4ir,  is  the  charge  density  on  the  wall,  and  is  the  potential 
difference  between  the  plane  of  <t/l  and  infinity,^as  shown  in  Fig.  1.  > 

If  e  «  e,  the  calculated  V  is  too  small.  Choosing  e<e  gives  a  larger  potential 
difference.  Obviously  e  can  be  adjusted  to  agree  with  experiment.  The  fact  that  a  single 
value  of  e  gives  fairly  good  agreement  with  experiment  for  a  wide  range  of  concentrations 
and  values  of  E  indicates  that  there  is  a  physical  basis  to  this  ad  hoc  procedure.  However,  as 


will  be  shown,  the  view  that  the  solvent  structure  need  be  taken  into  account  only  in  a  single 
layer  is  too  extreme. 


ION-DIPOLE  MODEL 

A  more  realistic  model  of  the  solvent  is  a  collection  of  hard  spheres  with  imbedded 
point  dipoles.  Although  lacking  as  a  model  of  water,  it  is  certainly  a  better  model  solvent 
than  a  dielectric  continuum  with  or  without  a  discontinuity  in  the  dielectric  constant  near  the 
wall.  Our  model  electrolyte  is  then  a  mixture  of  charged  hard  spheres  of  diameter  a  and 
charge  z^e  (where  Zj  is  the  valence  of  ionic  species  i  and  e  is  the  magnitude  of  the  electronic 
charge)  and  dipolar  hard  spheres  of  diameter  a  and  dipole  moment  /x.  We  could  consider  the 
ions  to  be  of  different  diameter.  We  will  not  do  this  firstly  because  there  has  been  little 


theoretical  work  on  the  theory  of  ions  of  unequal  size  and  secondly  because,  as  we  shall  see. 


the  explicit  solvent  term  dominates  and  so  the  effect  of  differences  in  ionic  diameters  is 


Let  p;«Nj/V  (V  is  the  volume)  be  density  of  the  ions  of  species  i  and  Pj^aN^/V  be 


the  density  of  the  dipoles.  Further,  let 


Because  of  charge  neutrality. 
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Near  the  wall  the  electrolyte  will  be  inhomogeneous.  Thus,  there  is  a  density  profile  for  the 
ions  which  we  may  write 

Pi(x)  »  Pjgjfx) 

»  Pj[hi(x)  +  1]  .  (3) 

where  x  is  the  distance  from  the  wall.  Similarly,  the  density  profile  of  the  dipoles  is 

p<,(x,Q)  «  pjg„(x,Q)  ,  (4) 

where  Q  specifies  the  orientation  of  the  dipoles.  Now, 

gj(x,fl)  »  g®(x)  +  /?  Ahj(x)  cos  d  +  ,  (5) 

where 

g^fx)  -  1  +  hS(x) 

"  47  /  -  (6) 

where  dfl  is  the  differential  element  of  solid  angle,  and 

Ahj(x)  «  ^  J  gd(x.n)  cos  0  dO  ,  (7) 

where  6  is  measured  from  a  normal  from  the  wall. 

If  —E/4v  is  the  charge  density  on  the  wall,  then  the  total  amount  of  charge 
accumulated  near  the  wail  must  be  equal  in  magnitude,  but  opposite  in  sign  to  this  charge 
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density.  Hence, 

®  S  «  E/4ir  . 


The  potential  difference  across  the  interface  is 


4ire2  ZjPi  f"  t  hj(t)dt 
i-1  0 


where  p  =  I  p  I . 


(8) 


THE  MEAN  SPHERICAL  APPROXIMATION 

Having  specified  the  model  and  defined  the  relevent  functions,  an  approximation  must 
be  introduced  to  obtain  results.  First,  we  use  the  mean  spherical  approximation  (MSA) 
where  the  core  condition  g(x)sO  inside  the  ionic  or  solvent  cores  is  combined  with  the 
approximation  that  outside  the  cores,  the  direct  correlation  function  is  equal  to  the  negative 
of  the  pair  potential  divided  by  kT  (k  is  Boltzmann's  constant  and  T  is  the  temperature). 

The  MSA  integral  equation  has  been  solved  for  the  ion-dipole-wall  system.^'^  The 
general  solution  is  implicit  and  has  not  yet  been  examined  fully.  However,  in  the  limit  of  low 
concentrations  of  ions,  the  MSA  results  are  explicit.  In  this  low  concentration  limit, 

ItM  -  god)  +  ^  , 


(10) 
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where  gQ(x)  is  the  density  profile  (normalized)  for  hard  spheres  near  a  hard  wallj  /3»l/kT, 
and 


41r^e^ 


2  zfpi 


(11) 


is  the  Debye  screening  length. 


Apart  from  the  term  g()(x),  which  makes  no  contribution  in  either  Eq.  (8)  or  (9) 
because  of  charge  neutrality,  Eq.  (10)  is  the  solution  of  the  MSA  for  a  system  of  charged 
hard  spheres  in  a  dielectric  medium  of  dielectric  constant  e.  In  other  words,  at  low 
concentrations  in  the  MSA,  the  ionic  terms  have  been  decoupled  from  the  dipolar  terms. 


In  the  MSA,  g®(x)  is  the  hard  sphere/hard  wall  profile,  go(x)-  In  the  limit  of  low 
concentrations,  Ah^(x)  satisfies 

>2 


Ahd(x)  -  SzipJj*"  hi(t)dt 

^  P‘3  i  (  *“*di 
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where 


(13) 


®di  2  ’ 

(a,b)  means  that  the  minimum  of  a  and  b  is  to  be  used,  and 

)33.2n  -  1  +  (  -  1)"  b2/(3-2")  . 
The  parameter  b2  (and  thus,  ^6*  related  to  e  by 

/>!  ■ 


(15) 


(14) 


The  MSA  values  for  the  two  integrals  which  contribute  to  V  in  Eq.  (9)  are 

4«2^,„/;ih,(t)d.  -  f[i  +  f]  (IS) 

and 

^  Pd/i  fjlb^(t)dt  -  -Lp.  [l  +  ■|((7-ad/\)]  f  .  (16) 

where  Thus, 

It  is  to  be  noted  that  there  has  been  a  cancellation  between  terms  of  order  1  and  (e-l)/e  to 
produce  a  term  of  order  1/e.  If 

•  0  , 
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we  obtain 


V  « 


(18) 


which  is  the  result  for  charged  hard  spheres  in  a  uniform  dielectric  medium.  Thus,  the 
contribution  due  to  the  discrete  nature  of  the  solvent  is  (e-1)  o^E/2Xe.  This  is  formally 
identical  to  Eq.  (1).  In  fact,  if 


e 


(19) 


Eqs.  (1)  and  (17)  are  identical  (if  the  MSA  is  used  to  calculate  <^q).  However,  the 
underlying  physical  pictures  are  quite  different.  In  Eq.  (1),  the  solvent  is  presumed  to  be 
affected  by  the  wail  only  within  an  inner  layer  of  thickness  a/2  whereas,  as  is  seen  from 
Fig.  2,  the  orientational  effects  of  the  dipoles  near  the  wall  extend  several  layers  into  the 
blind.  It  is  not  meaningful  to  talk  about  an  inner  layer  capacitor  in  series  with  a  diffuse  layer 
capacitor. 


Equation  (14)  can  be  rewritten  as 

16e  «  +  1)“*  .  (20) 

If  e  is  set  equal  to  the  bulk  dielectric  constant  of  water  (e~78),  Xa2.6S.  Thus, 
(e-l)/X-29,l.  A  reasonable  value  for  water  is  cTjj— 2.76A.  The  solvent  term  (e-l)oj/X 


dominates  so  that  the  precise  value  used  for  a  is  not  critical.  For  simplicity,  we  may  use 


s 


Results  obtained  form  Eq.  (17)  for  the  differential  capacitance, 


C  =  ii. 
^  ~  Ajt  dv 


e 

4ir 


■*f( 


a  + 


£-1 


(21) 


compared  with  experiment^  in  Table  1.  In  as  much  as  the  MSA  in  a  linearlized  theory,  it  is 
expected  to  be  useful  only  for  small  charge  on  the  electrode.  As  a  result,  the  comparison  in 
Table  1  is  confined  to  this  region.  The  agreement  is  good  considering  that  no  adjustable 
parameters  are  employed. 


A  SIMPLE  EXTENSION  OF  THE  MSA 

Well  known  extensions  of  the  MSA  exist.  Examples  are  the  hypemetted  chain  (HNC) 
approximation  and  more  simple  LHNC  and  QHNC  approximations.  The  LHNC 
approximation  is  a  linearized  theory  and  presumably  shares  many  of  the  deficiencies  of  the 
MSA.  The  others  are  more  complex  and  lead  to  equations  which  are  difficult  to  solve  even 
numerically. 


In  addition  to  the  numerical  difficulties,  there  e.\ists  another  problem.  We  have  seen 
that  in  the  MSA.  there  is  a  cancellation  between  terms  of  order  1  and  (e-l)/e  to  produce  a 
term  of  order  1  e  in  the  potential.  This  was  no  problem  in  the  MSA  since  all  the  terms  were 
analytic.  However,  in  a  theory  where  the  analogues  of  Eqs.  (15)  and  (16)  are  obtained 
numerically  and  where  e^SO,  a  one  percent  error  in  either  term  may  produce  a  hundred 
percent  error  in  V,  Clearly,  a  theory  must  be  formulated  so  that  this  cancellation  is  effected 
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before  the  commencement  of  numerical  calculations  if  a  reliable  result  for  V  is  to  be 
obtained. 

In  the  MSA,  Eq.  (12)  is  coupled  with  the  solution  of  the  MSA  for  hj(x)  for  charged 
hard  spheres  in  a  uniform  continuum  dielectric  medium.  A  simple  procedure  would  be  to 
retain  Eq.  (12)  but  to  use  some  more  general  expression  for  hj(x)  such  as  the  nonlinear  GC 
theory  or  the  HNC  theory  for  charged  hard  spheres  in  a  dielectric  background. 

In  general,  h{(x)  would  have  to  be  determined  numerically.  However,  as  long  as  the 
hj(x)  satisfy  charge  neutrality,  Eq.  (8),  the  potential  is  given  by 

V  -  ^  X  Vl4’  ^  .  (2J) 

The  desired  cancellation  between  terms  of  order  1  and  (e*l)/e  has  been  accomplished. 


Values  for 


^0  -  2  Vi/"  “  <^/2)  hi(t)dt ,  (21) 

®  i  •'*/2 

calculated  from  the  MSA,  GC  and  HNC  theories  for  charged  hard  spheres  in  a  dielectric 
medium,  are  compared  with  the  simulation  results  of  Torrie  and  Valleau^  in  Fig.  3.  The 
MSA  result  is  useful  only  for  small  charge  densities  on  the  wall.  The  GC  theory  results  are 
appreciably  larger  than  the  simulation  results.  On  the  other  hand,  the  HNC  results  are  in 
quite  good  agreement  with  the  simulation  results. 
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For  the  total  potential  V,  and  thus  the  differential  capacitance,  the  explicit  solvent 
contribution.  (e-l)E(7j[/2/\e,  dominates.  Hence,  as  is  seen  in  Fi^.  i,  the  differences  in  Cj 
between  the  use  of  the  GC  and  HNC  results  for  hj(x)  in  Eq.  (21)  are  smaller. 


SUM\L\RY 

Equation  (12)  provides  a  theoretical  basis  for  Eq.  (1).  The  potential  is  formally 
identical  that  obtained  from  a  collection  of  charged  hard  spheres  in  dielectric  medium  near  a 
charged  hard  wail  with  an  inner  layer  whose  properties  are  different  from  the  bulk. 

However,  the  physical  interpretation  is  quite  different.  The  effect  of  the  waU  on  the  solvent 
is  not  confined  to  an  inner  layer  of  thickness  0/2  but  extends  for  several  days  into  the  fluid. 
In  addition,  no  artificial  semi*macroscopic  concept,  such  as  a  local  dielectric  constant,  need 
be  introduced.  The  explicit  solvent  contribution  to  V  arises  from  microscopic  orientational 
effects  on  the  dipoles  near  the  wall. 

The  agreement  with  experiment  is  qualitative.  To  obtain  belter  agreement  with 
experiment  (e-l)/-\  would  have  to  be  treated  as  a  parameter.  The  charge  in  (e-l)/,\  would 
not  be  large  and  could  be  ascribed  to  polarization  or  some  other  phenomenon  neglected  in 
the  model. 

One  defect  of  Eq.  (22)  is  that  for  a  binary  electrolyte  for  which  I  Zj  I  -z.  the  resulting 
capacitance  is  synunetric  under  a  charge  in  the  sign  of  the  charge  on  the  electrode.  Allowing 
the  ions  to  be  of  different  size  does  not  help  much  since  the  explicit  solvent  term  dominates. 
Presumably,  the  asymmetry  comes  from  the  solvent.  Effects  such  as  a  nonspherical  solvent 
molecule,  specific  adsorption,  and  image  forces  should  be  uncorporated  into  the  model. 


..m 
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TABLE  1 


Differential  capacitance  (f/m^)  in  the  limit  of  zero  charge. 
(The  theoretical  values  are  for  (t--Ojj=s2.76A,  I  Zj  I  =z*l, 
T=298K,  and  e-78.4.  The  experimental  values  are  for  NaF.) 
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;  a  function  of  llic  charge  on  the  electrode  for 
6=78.4.  The  solid  and  broken  curves  give  the  results 
theories,  respectively,  fur  h|(x). 


